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ABSTRACT 



We have analyzed the bivariate distribution of galaxies as a function of 
ultraviolet-optical colors and absolute magnitudes in the local universe. The 
sample consists of galaxies with redshifts and optical photometry from the Sloan 
Digital Sky Survey (SDSS) main galaxy sample matched with detections in the 
j> ■ near-ultraviolet (NUV) and far-ultraviolet (FUV) bands in the Medium Imaging 

Survey being carried out by the Galaxy Evolution Explorer ( GALEX) satellite. 
Q\ • In the (NUV — r) 01 vs. M r0 i galaxy color-magnitude diagram, the galaxies 
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separate into two well-defined blue and red sequences. The (NUV — r)o.i color 
distribution at each M r .i is not well fit by the sum of two Gaussians due to an 
excess of galaxies in between the two sequences. The peaks of both sequences 
become redder with increasing luminosity with a distinct blue peak visible up to 
M r fi.i ~ —23. The ro.i-band luminosity functions vary systematically with color, 
with the faint end slope and characteristic luminosity gradually increasing with 
color. After correcting for attenuation due to dust, we find that approximately 
one quarter of the color variation along the blue sequence is due to dust with the 
remainder due to star formation history and metallicity. Finally, we present the 
distribution of galaxies as a function of specific star formation rate and stellar 
mass. The specific star formation rates imply that galaxies along the blue se- 
quence progress from low mass galaxies with star formation rates that increase 
somewhat with time to more massive galaxies with a more or less constant star 
formation rate. Above a stellar mass of ~ 10 10 ' 5 M , galaxies with low ratios of 
current to past averaged star formation rate begin to dominate. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters - 
galaxies: luminosity function — galaxies: statistics — galaxies: surveys — ul- 
traviolet: galaxies 



1. Introduction 



Galaxies exhibit bimodal distributions in a number of observed properties. The bimodal- 
ity in g a laxy morphologies formed the basis of the original galaxy classification scheme of 
Hubble! ( 119261 ). The c olors and luminosities of galaxies have been lo ng known to correlate 



with morphology (e.g. Ide Vaucouleurslll96ll ; IChester fc Roberts! 1 19641 ) with ellipticals being 



predominantly red and spirals and irregulars blue. 

More recently, large statistical samples of galaxies have become available, allowing us 
to investigate the bimodality of galaxies in a much more quantitative way. In particular, 
the bimodality appears quite strongly in the galaxy (u — r) color distribut ion which con- 
sists of two peaks with a minimum in between them at {u — r) ~ 2.1 — 2.2 (IStrateva et al. 



200ll ). Galaxies in the red peak tend to be predominantly morphologically early- type and 
high surface brightness galaxies while those in the blue peak are dominated by morpho 



logically late-type galaxies with lower surfa ce brightness ( IStrateva et al.l l200ll ; iDriver et al. 



2006; 


Blanton et al. 


2003b: 


Ball et al. 


2006) 



vs. 



M r color- magnitude diagram (CMD). The galaxies in their sample separate into blue 
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and red sequences with the distribution in color at each absolute magnitude well-fit by the 
sum of two Gaussians. The mean color as a function of M r for each sequence consists of an 
overall reddening with increasing luminosity with a steeper transition in the average color 
and width of both sequences at a stellar mass of ~ 2 x 10 10 M . 

In addition to mass, one of the most important other factors suspected of contributing 
to the galaxy bimodality is the environment. While it has long been kno wn that the mor- 
phologies of galaxies are correlated with the local density (jDresslerlll980l ). the dependence 
of galaxy colors and luminosities with local density is complicated. Although the ratio of 
the number of red to blue galaxies varies strongly with the local den sity, the mean colo r 
of the blue and red sequences varies relatively little with environment (IBalogh et al.l 12004 ). 
On the other hand, the luminosity of blue sequence galaxies is nearly independent of envi- 
ronment while both luminous and faint red galaxies are found on average in higher density 
environments than intermediate luminosity red galaxies (IHogg et al.ll2003l ). 



The galaxy bimodality has also begun to be investigated based upo n large samples of 



galaxy spectra from the Sloan Digital Sky Survey (SDSS). In particular, iKauffmann et al. 



( 2003al ) developed a method that uses the Balmer absorption line index H<5a and the 4000 A 
break strength D n (4000) measured from the SDSS fiber spectra in the central 3" of each 
galaxy to constrain the star formation histories, d ust attenuation, a n d stell ar masses for 
their sample. Based upon these derived parameters, IKauffmann et al.l (l2003bl ) showed that 
galaxies tend to divide into two distinct groups around a ste llar mass of 3 x 1 10 M , 
similar to the transition mass noted in the optical galaxy CMD (IBaldry et al.ll2004l ). While 
galaxies below this mass tend to ha ve younger stellar popula tions, more massive galaxies 
tend to be older. In related work, iBrinchmann et al.l ( 2004 ) used the emission lines in 
the SDSS spectra to determine star formation rates (SFRs) for a large sample of SDSS 
galaxies. Using the specific star form ation rate, i.e. the current SFR divided by the stellar 
mass M*, IBrinchmann et al.l (12004 ) found that galaxies with 10 8 < M < 10 9 M have 
log (SFR/M*) = —9.6 to —10, values consistent with an approximately constant SFR with 
time. Above 10 10 M Q , the specific SFRs decline with mass, implying star formation histories 
increasingly weighted to much older ages. 

The evolution of the galaxy color-magnitude diagram out to z ~ 1 has begun to be 
explored (jWillmer et al.l 120061 ; iFaber et al.l 120061 ; iBlantonl 120061 ) . These results show that 
the galaxy bimodality is already in place at z ~ 1. However, the color of both sequences 



tends to become somewhat bluer with increasing redsh ift while the 



both red and blue galaxies shifts to higher luminosities dBlanton 
Based upon combining the DEEP2 and COMBO 17 surveys, 



2006 



uminosity function o f 



Willmer et al.l 120061) . 



Faber et al. 



hood ) argued that 



the number density of blue galaxies is more or less constant from z ~ 1 to z ~ 0, while the 
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number density of red galaxies has been increasing. iFaber et al.l (120061 ) proposed a scenario 
to explain their data in which some blue galaxies migrate to the red sequence as a result of 
gas-rich mergers that use up the remaining gas in an interaction-induced starburst. These 
galaxies then migrate up the red sequence by a series of gas-free mergers. 

The origin of the galaxy bimodality and corresponding transition mass of a fewxlO 10 
M is beginning to be understood theoretically. Based upon a semi-analytic model utilizing 
some simple prescriptions for gas cooling, star formation, and su pernova feedback coupled 
with mergin g histories of d ark m atter haloes, iMenci et al.l (120051 ) modeled the (u — r) vs. 
M r CMD of iBaldry et al.l (120041 ) . In their model, feedback from supernovae is ineffective 
at regulating star formation for galaxies above a certain threshold halo mass. In these 
massive galaxies all of the gas is consumed relatively quickly and results in a red sequence 
galaxy at zero redshift. Blue sequence galaxies, on the other hand, tend to come from less 
massive progenitors where supernovae feedback is effective at regulating star formation, thus 
allowing star formation to continue down to the present. While their model is successful 
at reproducing most of the optical CMD, it predicts too many blue galaxies at M r = —22 
compared to the observations. 



A different explanation for the origin of bimodality has been suggested by lDekel &: Birnboim 

10 12 M n , a shock is 



( 120061 ). According to this model, above a critical halo mass M. 



shock 



generated in the gas accreting onto the dark matter halo which heats most of the gas and 
prevents it from cooling and forming stars. In these massive haloes, star formation does 
happen at z > 2 due to cold gas that is able to penetrate the hot gas, leading to a burst of 
star formation, while for z < 2 heating from Active Galactic Nuclei (AGN) prevents gas from 
forming any more stars. This naturally leads to the most massive galaxies lying on the red 
sequence at z ~ 0. For galaxies residing in halos with masses less than 10 12 M , the gas is not 
shock heated, allowing cold flows to fuel star formation that is then regulated by supernova 
feedback. As a result, lower mass galaxies lie on the blue sequence and the location of the 
bright tip of the blue sequence is due to the onset of the shock in the accreting gas for more 
massive halos and the feedback from AGN. In this scenario, galaxies tend to move up the 
blue sequence with time until their masses go above M s h oc k, or they merge into another more 
massive halo with mass above M s h oc k, after which the gas in t he g alaxy is no longer al lowed 
to cool and star formation ceases. Both lCattaneo et al.l (120061 ) and lCroton et al.l (120061 ) have 
coupled semi-analytic models including the transition from shock heating to cold flows and 
feedback from supernova as well as AGN with the merging histories of dark matter haloes 
from N-body simulations. While the details of the modeling of the baryonic physics differs 
somewhat, both groups were able to reproduce the local galaxy CMD by tuning the various 
parameters affecting star formation and feedback in their models. 
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In this paper, we investigate the galaxy bimodality as revealed in the UV minus optical 
colors of a large sample of galaxies observed by both the Galaxy Evolution Explorer (GALEX) 
and the SDSS. While significant contributi ons to the UV luminosity can c ome from older 
evolved stars in red sequence galaxies (e.g. lYi et al.ll2005l ; iRich et al.ll2005l ). in general the 
UV light in galaxies is dominated by massive stars with main sequence lifetimes up to ~ 10 8 
yrs. As a result, the emerging UV luminosity is proportional to the recent star formation 
rate once corrected for light absorbed by dust (IKennicuttl Il998l ). The greater sensitivity 
of the GALEX bands to the recent star formation rate as compared to the SDSS u band 
would lead us to expect a greater separation between the red and blue sequences. While the 
measurements from the SDSS spectra are sensitive diagnostics of the stellar populations in 
galaxies, they are measured only in the central 3" of each galaxy, making somewhat uncertain 
aperture corrections necessary to account for the portion of each galaxy not sampled. While 
the UV data presented here is much more susceptible to dust attenuation than in the optical 
SDSS data, the UV measurements sample the entire galaxy and thus complement the SDSS 
measurements. 



2. Data and Analysis 



2.1. GALEX Data 



Th e UV data present ed here are derived from the GALEX Medium Imaging Survey 
(MIS) ( IMartin et al.l 120051 ). GALEX is a 50cm diameter UV telescope that images the 
sky simultaneously in both a FUV and a NUV band, centered at 1540 A and 2300 A, 
respectively. The field-of-view of GALEX is approximately circular with a diameter of 1?2 
and resolution of about 5'/5 FWHM in the NUV. The MIS pointings are chosen to overlap 
areas of sky with imaging and spectroscopy from the SDSS and consist of exposures of at 
least one to a few orbits with the mode of the exposure time distribution being 1700 sec. 
The dataset used in our analysis is taken from the union of the GALEX first data release 
(GR1) with the GALEX internal release IR1.1, a subset of which has been included in the 
second data release (GR2) publicly available from the GALEX archivefl The IR1.1 data was 
processed with a pipeline very similar to that used in the GR1 data and employed the same 
calibration as used in that release. D etails of the GALEX detect ors, pipeline, calibration 
and source extraction can be found in iMorrissey et al.l (120051 . 120071 ) . 



The GALEX pipeline uses the SExtractor program (IBertin fc Arnoutd Il996l ) to de 



1 The GALEX archive can be accessed from 



http://galex.stsci.edu/GR2/ 



- 6- 



tect and make measurements of sources in the images. Throughout this paper, we use 
the "MAG_AUTO" measurements output by SExtractor. These magnitudes are measured 
within elliptical apertures with semi-major axis s c aled t o 2.5 times the first moment of each 
source's radial profile, as first suggested by iKronl (ll980l ). 



Due to an error in the way in which the SExtractor parameters were set in the GALEX 
pipeline, the photometric errors for most of the sources, as reported in the GR1 and IR1.1 
catalogs, are underestimated. For each source we have calculated a more accurate statistical 
error in the total magnitude from the size of the MAG_AUTO aperture and the flat field 
response, exposure time, and sky background at the source position. In addition to the sta- 
tistical errors, we have added in quadrat ure an assumed zero-p oint plus flat field uncertainty 
of 2% in the NUV and 5% in the FUV (IMorrissey et al.ll2007l ). The errors increase from the 
zeropoint uncertainty at the bright end up to ~ 0.2 — 0.3 mag at 23rd mag in both bands. 



2.2. GALEX-SBSS Matched Sample 

The GALEX MIS catalogs were matched with the SDSS MPA/JHU DR4 value-added 
catalogs^ These catalogs consist of line and index measurements from the SDSS spectra as 
well as many derived quantities and are described in more detail in a series of papers on the 
star formation rat es, star formation histor ies , stellar masses, and metallicities of galaxies i n 
the local universe fcauffmann et aI.ll2003al fbL I Brinchmann et al.ll2004l ; iTremonti et al.ll2004l ) . 



For each GALEX pointing, SDSS sources within 0?6 of the GALEX field center were matched 
with the nearest GALEX source within a radius of 4". When concatenating together the 
catalogs for all the fields, we removed duplicate GALEX detections in the overlap regions 
between adjacent pointings by using the SDSS identification numbers (Plate ID, MJD, Fiber 
ID) and selecting the GALEX match closest to its field center. 

After matching the GALEX and SDSS data, we further restricted the sample with 
various cuts intended to generate a complete statistical sample which are summarized in 
Table [TJ For the SDSS photometry, we selected galaxies targetted for spectroscopy in the 
SDSS main galaxy sample with r-band magnitudes in the range 14.5 < r < 17.6. While th e 



nominal magnitude limit of the SDSS main galaxy sample is r = 17.77 (IStrauss et al.ll2002l ). 
in practice the limit varies as a function of position on the sky. After examining the galaxy 
number counts, we set the faint limit to r = 17.6 because the counts begin to turn-over below 
this level. While the median photometric error for the SDSS galaxies is only 0.03 mag, there 
are a small fraction with much larger errors. In order to remove these objects, we restricted 



2 The MPA/JHU value-added catalogs were downloaded from 



http://www.mps-garching.mpg.de/SDSS/ 
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the sample to galaxies with errors a r < 0.2 mag. In addition to the photometry, we further 
restricted the sample to those galaxies with redshifts z in the range 0.01 < z < 0.25 and 
with redshift confidence z con f > 0.67. 

In addition to the cuts on the SDSS data, we applied several cuts based upon the 
UV measurements. Since GALEX photometry and astrometry degrade near the edge of 
the detectors, we only included objects in the sample if their distance from the GALEX 
field center fovjradius is less than 0?55. GALEX detections were also required to have 
nuv-artifact < 1. This excludes from our sample galaxies that lie within regions expected 
to be contaminated by reflections from bright stars within the field. Areas of the images 
with nuv jarti j 'act = 1 are designed to flag regions where scattered light is predicted from 
bright stars just outside the field-of-view. Currently, this flag is set very conservatively and 
the vast majority of sources with this flag set are in regions that are free from scattered light 
issues. We therefore elected to ignore this flag for our sample. We also only included fields 
with exposure times greater than 1000 sec. The resulting overlap area between GALEX and 
the SDSS including all of the above cuts is 485.321 deg 2 in the NUV and 411.266 deg 2 in 
the FUV. The FUV sample has a somewhat smaller area because some of the fields have 
NUV data on l y. Th e procedure used to calculate the overlap area is similar to that used by 



Bianchi et al.l (120071 ). where a more detailed description can be found. 



For both the FUV and NUV samples, we included galaxies with apparent magnitudes 
in the range 16 < N(F)UV < 23. Based upon a series of artificial source tests, we have 
estimated that the data is greater than 80% complete for NUV, FUV < 23 mag. After 
applying the Galactic extinction correction, the UV magnitude limit varies across the sky. 
We account for this variation below when computing the volume densities. After applying 
all of the above cuts to the data, the NUV and FUV samples contain 26,281 and 18,091 
galaxies, respectively. The redshift distributions for the NUV and FUV samples are shown 
in Figures [T] and [21 In the figures, the solid black lines are the distributions of SDSS galaxies 
lying within the area observed by GALEX and satisfying the SDSS selection criteria while the 
dashed red lines show those galaxies which have a GALEX match falling within the GALEX 
selection criteria as well. The fraction of SDSS main sample galaxies with GALEX matches 
tends to decline with redshift, mainly due to increasing numbers of red galaxies falling below 
the GALEX detection limit. In Figure [3] the r-band magnitude is plotted as a function of 
color for the NUV and FUV samples. The dashed lines in the figure indicate the r-band and 
UV magnitude limits. For the NUV sample, the lack of galaxies with (NUV — r) > 6.5 at 
bright r magnitudes would argue that the GALEX MIS data is deep enough to fully sample 
the entire color distribution of galaxies in the local universe. Thus, the edge of the color 
distribution is likely real and not a selection effect. On the other hand, for the FUV sample 
in the right hand panel of Figure [31 the data go right up to the selection limit at the red end. 
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In the FUV sample, the red edge of the color distribution thus reflects the UV magnitude 
limit. 

The fraction of SDSS main sample galaxies that lie within 0?55 of a GALEX MIS field 
center and that have a GALEX match within our 4" search radius is shown in Figure H] as 
a function of r magnitude. In the NUV sample, the completeness is roughly constant at 
about 90% down to r ~ 16.5. Fainter than this magnitude, the completeness begins to drop 
off. There are very few galaxies with colors redder than (NUV — r) « 6.5. At the GALEX 
magnitude limit of 23 mag, a galaxy with this red a color would have r = 16.5. Thus, the 
fraction of SDSS galaxies with a GALEX match begins to drop fainter than this r magnitude 
due to increasing numbers of red galaxies falling below the GALEX detection limit. For the 
FUV sample, the reddest galaxies have (FUV — r) « 7.5 which corresponds to r = 15.5 at 
the limiting magnitude of FUV = 23. As expected, the fraction of SDSS galaxies with an 
FUV match begins to decline at about this r magnitude. 

In both the FUV and NUV samples, the match completeness does not reach 100% 
at bright r magnitudes. We have visually inspected all SDSS main sample galaxies with 
14.5 < r < 15.5 that were observed by GALEX. For almost all of these galaxies, there is 
a galaxy visible in the GALEX images. However, the UV center measured by the GALEX 
pipeline for these non-matches lies more than 4" from the SDSS position. Sometimes the 
GALEX pipeline breaks the galaxy into more than one fragment, none of which are coincident 
with the SDSS position. In other cases, especially if the galaxy has a low UV surface 
brightness, the center can be offset from the SDSS position by more than our search radius 
even if the GALEX pipeline detects the galaxy as a single object. We assume that the level 
portion of the match completeness curves in Figure H] gives the intrinsic completeness for 
GALEX detections of the SDSS main galaxy sample. The values we adopt are 0.91 and 0.80 
for the NUV and FUV, respectively. 

While the completeness of the SDSS photometric sample is nearly 100%, some fraction 
of galaxies that meet the SDSS main galaxy sample selection do not have a redshift measured 



(IStrauss et al.ll2002l ). Although some galaxies do not have a redshift due to low signal-to- 
noise in their spectra, the majority of targeted galaxies without redshifts are missed due 
to the constraint that SDSS spectroscopy fibers can not be placed closer than 55" to one 
another. Some of these missed galaxies can be observed in neighboring plates if that region 
of sky is covered by more than one plate. While the exact completeness is determined by the 
precise geometry of the spectroscopic plates, the result is that the spectro scopic completeness 



of the SDSS main galaxy sample is 92-94% for the early release data (IStrauss et al.l 12002 



Blanton et al.ll200ll ). We have adopted a spectroscopic completeness of 0.9. Multiplying the 



GALEX-SDSS match completeness by the SDSS spectroscopic completeness, we estimated 
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the total completeness of our sample to be 0.82 in the NUV and 0.72 in the FUV. In 
calculating the volume densities below, we correct the number counts by these factors, (i.e. 
the factor / in equations (jSJ) and (J5J) below). 

In addition to the spectroscopic incompleteness, there is an additional surface brightness 
selection that is imposed on the SDSS main ga laxy sample. A s a pa rt of their study of the 



luminosity function of low luminosity galaxies, iBlanton et al.l (120051 ) investigated the com- 
pleteness as a function of surface brightness and found that the SDSS spectroscopic galaxy 
sample is greater than 90% complete above a half light surface brightness of /i 50ir = 22.4 
mag arcsec -2 with the completeness dropping to 50% at /i5o, r = 23.4 mag arcsec -2 . Since lu- 
minosity and surface brightness are correlated, the surface brightness s election preferentially 



selects against dwarf galaxies. For galaxies brighter than M r = —18, IBlanton et al.l (120051 ) 
have a fit a gaussian to the surface brightness distribution in a series of absolute magnitude 
bins and have used this model to extrapolate the number of galaxies likely missed due to the 
surface brightness selection at fainter luminosities. The fraction of galaxies missing from the 
sample increases from near zero at M r = —18 to approximately 40% at M r = —16. However, 
these low luminosity, low surface brightness galaxies do not make a significant contribution 
to the total luminosity density. Even after correcting for the surface brightness incomplete- 
ness, about 90% of the r-band luminosity density is due to galaxies with M r < —17. Thus, 
the surface brightness selection is unlikely to affect our results for bright galaxies. However, 
we may be underestimating the number density of galaxies with —18 < M r < —16 by up to 
40%. 



2.3. Absolute magnitudes and volume densities 

We computed absolute magnitudes for our sample galaxies using, for example for the 
r-band, 

M rA1 = m r -5\ogD L - 25 - K^ r (z) + (z - 0.1)Q (1) 

where M r .i is the absolute magnitude, m r is the extinction corrected r-band magnitude, 
Dl is the luminosity distance in Mpc, K i r (z) is the K-correction needed to account for 
the shifting of the galaxy SEDs with respect to the filter bandpass, and Q is a term to 
account for luminosity evolution in units of magnitudes per redshift. A positive value for 
Q means that galaxies get brighter with increasing redshift. Similar equations were used 
for the other bands. For calculating the luminosity distance we assumed a Hubble constant 
Hq = 70 km s _1 Mpc -1 and a flat universe with matter density relative to the critical density 
of Q m = 0.3 and dark energy density of = 0.7. We calculated the K-correcti o ns usin g the 



K_CORRECT program, version 4.1.4, originally developed by IBlanton et al.l (j2003al ) and 
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now extended to handle GALEX data (IBlanton fc Roweia 120071 ). Given the redshift of a 
galaxy, the K_CORRECT program finds the linear combination of a set of template galaxy 
spectra that best reproduces the observed colors of a galax y. The templates have already 
been determined as described by IBlanton fc Roweisl (120071 ). The coefficients in the linear 
combination for each galaxy were used later when determining the maximum volume out 
to which a given gala xy could be observed and be detected in the sample. As suggested by 
Blanton et al.1 (j2003b| jd). we minimized the errors due to the K-corrections by correcting all 
of the galaxies to bandpasses shifted by z = 0.1, a redshift near the median value for our 
sample. We denote absolute magnitudes and colors in this system with the subscript 0.1. 

Probably the most uncertain term in equation ([1]) is the evol ution term Q. When 



calculating the optical galaxy luminosity functions from SDSS data, IBlanton et al.1 (l2003bl ) 
fit for both luminosity and density evolution even within the relatively small redshift range 
sampled by the SDSS data. While the results for number and density evolution are highly 
correlated, the best fit luminosity functions are consistent with no nu mber density evolution 
but significant luminosity evolution. In particular for the ro.i band, IBlanton et al.1 (l2003bl ) 
found Q = 1.62 ± 0.3, corresponding to a diff erence of 0.4 mag over the redshift range in 
our sample of 0.01 < z < 0.25. As shown by IBlanton et al.1 (j2003bl ). neglecting evolution 
can lead to significant distortions in the shape of the luminosity function. We have adopted 
Q — 1.6 for all bands considered here. This means that we have implicitly assumed that 
galaxies evolve in luminosity but not in color or in number density. 

A value for Q o f 1.6 is roughly consistent with other determinations in the UV as well. 
Schiminovich et al.1 (120051 ) investigated the evolution of the UV luminosity density Pfuv with 
redshift and found for z < 1 that Pfuv oc (1 + z) 2 - 5±0 - 7 , corresponding to Q = 1.9 ± 0.5, 
assuming that the increase in ppuv is due entirely to luminosity evolution. Based upon a 
sam ple of GALEX galax i es matched w i th red shifts from the Two Degree Field Redshift Sur- 
vey (IColless et al.ll200ll ). iTreyer et al.1 (120051 ) measured evolution of the luminosity function 
characteristic magnitude M* of AM^ UV « 0.3 ± 0.1 mag between z = 0.05 and z = 0.15, 
equivalent to Q = 3 ± 1. These other measurements would tend to favor a somewhat larger 
evolution in the UV with redshift as compared to the r-band. Indeed, galaxies most likely 
become bluer with increasing redshift as a result of the average star formation rates of galax- 
ies increasing with redshift. For example, comparing a low-z sample from the SDSS and a 



sample of galaxies at z ~ 1. IBlanton! (120061 ) found that the blue sequence becomes bluer by 
0.3 mag in (u — r) while the red sequence beco mes bluer by only 0.1 mag. In another study 



of the evolution of the blue and red sequences, iFaber et al.1 (120061 ) analyzed the evolution of 
galaxies in the (U — B) vs. Mg CMD out to z ~ 1. They found that Mg becomes brighter 
by ~ 1.3 mag out z ~ 1 for both blue and red sequence galaxies. On the other hand, the 
luminosity function normalization <ft* for blue galaxies was found to be roughly constant 
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while that of red galaxies is increasing. 

In order to asses the effect of color evolution on our results, we have recomputed the 
absolute magnitudes and volume densities described below for the NUV sample assuming 
evolution in the r-band of Q r = 1.6 and in the NUV of Qnuv — 3. These choices correspond 
to a decrease in the (NUV—r) color of a galaxy across our redshift range from 0.01 < z < 0.25 
of 0.34 mag. We have compared many of the results presented in the following sections with 
and without allowing for color evolution. The morphology of the color-magnitude diagram 
remains largely unchanged with the peaks and widths of the red and blue sequences nearly the 
same. The largest effect of including color evolution is on the volume density of luminous 
blue galaxies since these galaxies are detectable across the entire redshift range and thus 
would have the largest color correction. Including color evolution also tends to increase 
somewhat the volume density of galaxies with very blue colors (NUV — r) < 1. Specifically, 
including color evolution would increase the luminosity density of the bluest galaxies with 
< (NUV — r) < 1 by 0.2 dex in the r-band and by 0.3 dex in the NUV. Due to the large 
uncertainties remaining in the evolution of galaxy colors with redshift and the relatively 
minor effect it has on our results, we decided to neglect color evolution in our analysis. 



When correcting for Galactic extinction, we assumed the lCardelli et al.l (119891 ) extinction 
law with R v = A V /E(B - V) = 3.1. For the SDSS bands, the ratio of A(X)/E(B - V) is 
5.155, 3.793, and 2.751 for u, g, and r, respectively, while for the FUV the ratio is 8.24. 
Due to the presence of the 2175 A bump in the Galactic extinction law, the extinction in 
the NUV band is no longer strictly proportional to the reddening E(B — V). In order 
to quantify the effect this ha s on o ur extinction corrections, we used a small set of 42 
SEDs from lBruzual fc Charlotl (120031 ) that span a representative range of galaxy SEDs from 



quiescent ellip t icals t o rapidly star-forming galaxiesfj For each intrinsic SED, we applied the 



a repi 
ies@] 



Cardelli et al.l (119891 ) extinction law and then computed the resulting NUV AB magnitude 



as a function of E(B — V). For each SED, we fit a a quadratic function of E(B — V): 

A NUV = a x E(B -V) + a 2 E(B - V) 2 . (2) 

For galaxies with some recent star formation, a\ = 8.24 and a 2 = —0.67 while for older 
galaxies with little or no recent star formation and a SED that falls steeply in the UV with 
decreasing wavelength, ai is slightly smaller and lies in the range 7.5 — 8.0. For 97% of our 
sample E(B — V) < 0.1, and thus the quadratic term can be safely neglected. In addition, 
the maximum difference in the adopted value for A(NUV) for the range of values for ai 
among the 42 SEDs, is only 0.07 mag at E(B — V) =0.1. Therefore, we assume the value 
A(NUV)/E(B -V) = 8.2 for all of our calculations. 



3 The 42 SEDs can be found at |http://www.lam.oamp.fr/arnouts/LE_PHARE .html 
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As argued above, the fraction of SDSS main sample galaxies with GALEX detections 
is a strong function of the galaxy color. In Figures [5] and [61 we plot contours of the fraction 
of SDSS galaxies with a GALEX match for the NUV and FUV samples, respectively. In 
both samples, the completeness for blue galaxies is more than 90% while the completeness 
begins to drop for galaxies with (g — r)o.i > 0.8. For the NUV sample the completeness 
along most of the red sequence is in the range 30 — 60% while the completeness of the red 
sequence in the FUV sample is lower and in the range 10 — 40%. It is important to note 
that this drop-off in the fraction of galaxies with a GALEX match is due to the GALEX 
magnitude limit and was taken into account below when computing the volume densities of 
galaxies as a function of absolute magnitude and color. 



We used the V max method (jSchmidtlll968l ) to determine the volume densities of galaxies 



in our samples. The value of V max for each galaxy is given by the maximum volume within 
which the galaxy could have been included in the sample, given our selection limits listed in 
Table [TJ We computed a separate V max for the FUV and NUV samples. First, we computed 
Kq.i(z) for 0.01 < z < 0.25 for each galaxy using the best-fit SED derived from the output of 
the K_CORRECT program. Next, we used equation (OQ) to define a maximum and minimum 
redshift for each galaxy in each band by replacing the apparent magnitude in equation ([TJ 
by the magnitude limits listed in Table [H Then we computed a combined minimum and 
maximum redshift using 

Zmax Wl'il7'(.Zr,rnax j ^UV,max j 0.25) , (3a) 
Zmin TTlClX{Z r)m i n ) ZUV,mim 0.01). (3b) 

Since we have assumed a cosmology with no overall curvature, we calculated the volume 
between z min and z max for each galaxy as 

v - A ( 71 Y ( Dh ^ Zmax ^ D dz min f \ 

^ ~ 3 \l80J V(l + Zmaxf (1 + Z zmm f) ' 1 J 

where A is the solid angle in deg 2 from which the sample was drawn. The V max values can 
then be used to generate the number densities of galaxies as a function of any variables. 
For example, in order to generate a volume-corrected galaxy CMD for the NUV sample, 
we computed the number density of galaxies as a function of absolute magnitude and color 
using 

* W , o ,,(W-r) o .0 = ^£^, (5) 

where $ gives the number density of galaxies in units of Mpc -3 mag -2 , AM is the width of 
each bin in absolute magnitude, AC is the width of each bin in color, / is the inverse of the 
sample completeness estimated above in §2.2, and the sum is taken over all galaxies within 
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that particular color and absolute magnitude bin centered on M rj o.i and (NUV — r)o.i- The 
corresponding uncertainty in each bin due to counting statistics is given by 

^, al ,(w- r )o.o = ^(£^) 1/2 . (6) 



3. Results 

3.1. The Galaxy Color Magnitude Diagram 

The number of galaxies in our (NUV — r) .i and (FUV — r) 01 CMDs are plotted in 
Figures [7] and El respectively. In both plots, the data are plotted as contours where the 
density of points is high while the locations of individual galaxies are plotted where the 
density is low. The uncertainty in the colors in both samples is dominated by the errors 
in the UV measurements. Thus the errors as a function of position in the CMDs are most 
strongly correlated with color. The median error as a function of color is plotted in Figures [7] 
and M along the left-hand side of each figure. For blue galaxies, the uncertainty is dominated 
by the zero-point uncertainty in the GALEX data. For red galaxies, the errors span a larger 
range from 0.1 — 0.4 mag with a median of about 0.2 mag. 

The corresponding volume densities of galaxies as a function of position in the CMD 
are plotted in Figures M and (TDl where the weighting was derived from the V max values as in 
equation (jSj) with AM = 0.5 mag and AC = 0.2 mag. The peak of each sequence from the 
Gaussian fits described below are over-plotted as the dashed lines in the (NUV — r) diagram. 
The volume densities and the errors for the (NUV — r) diagram are given in Tables [2] and 
[3] while those for the (FUV — r) diagram are given in Tables H] and 

As we have argued in §2.2, the red edge of the color distribution in the NUV sample 
is real whereas the red edge of the FUV sample is a selection effect due to the FUV flux 
limit. This is reflected in the morphology of the galaxy distributions in Figures [9] and [TOj 
where the distribution turns over for the reddest colors in the NUV diagram and does not 
turn over entirely in the FUV diagram. Throughout the remainder of this paper, we focus 
on the NUV diagram. 

In both the FUV and NUV diagrams, the galaxies separate into two well-defined se- 
quences in addition to a po pulation of galaxie s that lie in between. As has been noted 



before in optical CMDs (e.g. iBaldry et al.l 120041 ). the most luminous galaxies are on the red 



sequence, while both sequences become redder with increasing luminosity. In contrast to the 
optical (u — r) CMD, the blue sequence does not appear to merge at the bright end with the 
red sequence. 
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An alternative view of the NUV sample is shown in Figure [IT] where the volume density 
of galaxies is plotted as a function of the NUV luminosity. The sample reaches significantly 
fainter NUV absolute magnitudes for the red galaxies due to the SDSS r-band selection. 
Thus, the slope for the faintest NUV absolute magnitudes included in our sample as a 
function of color is a selection effect. As in Figure [9j the sample separates into blue and 
red sequences. However, there is little, if any, trend of color with Mnuv,o.i along either 
sequence. This is consistent with the conclusions from studies in the optical which indicate 
that one of the most important factors determining the evolution of a galaxy is its mass, 
which is much more closely related to the r-band luminosity than to the NUV luminosity 
( IKauffmann et al.l l2003bl ; iBrinchmann et al.l 120041 ) . 



3.2. Color Distributions as a Function of M. 



r-,0.1 



The volume-corrected number density of galaxies as a function of (NUV — r) .i color 
is plotted in Figures [TJHH] in 0.5 magnitude wide bins of M rj0 .i- The error bars are the 
statistical errors only calculated using equation ([6]). Except for the most luminous bin, there 
is both a red and blue peak visible in each panel. Similar to previous optical CMDs, the red 
sequence dominates in the brighter bins. The red and blue sequences reach approximately 
equal strengths around M rj0 .i = —21.75 with the blue sequence becoming dominant at fainter 
luminosities. The relative number of red sequence galaxie s reaches 50% at a bout the same 
luminosity when dividing galaxies using the (u — r) color (IBaldry et al.ll2004j ). 



Following iBaldry et al.l (120041 ). we attempted to fit Gaussians to the red and blue peaks 
in the color distributions in each M r0 .i bin although we employed a somewhat different 
methodology. We fit a single Gaussian of the form (1/ '\f2na 2 ) exp {—((NUV — r) .i — fi) 2 /2a 2 } 
separately to the red and blue seque nces. In order t o select points on the red and blue se- 
quences, we have utilized the fit from lYi et al.l (120051 ) to the (NUV — r) color as a function 
of M r for a sample of morphologically selected early-type galaxies: (NUV — r) = f(M r ) = 
1.73 — 0.17M r . For the purposes of fitting a Gaussian to each sequence, we defined the 
red sequence as the points with (NUV — r) .i > f(M r0 i) — 0.5 and fit a Gaussian to 
these points using the IDL routine GAUSSFIT, a program that computes a non-linear least 
squares fit to the data. Similarly, we fit a Gaussian to the blue sequence for points with 
(NUV — r) .i < /(M rj0 .i) — 2.0. These color limits are plotted as the dashed red and blue 
lines in each of the color distributions shown in Figures IT2T-IT41 whereas the sum of the best- 
fitting Gaussians is plotted as the solid line. In contrast to the optical (u — r) CMD, a 
double-Gaussian does not provide a good fit to the data. Although the Gaussians provide 
a reasonable fit to the blue edge of the blue sequence and the red edge of the red sequence, 
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their sum falls well below the data in the region between the two sequences. 



For the star-forming galaxies in the blue sequence, it is difficult to generate galaxies with 
extr emely blue colors, i .e. (NUV — r) 01 < 1, except with large starbursts or very young ages 
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Treyer et al.lll998l ). Clearly, such objects are rare in the local universe. The skew of the 
blue sequence in the (NUV — r)o.i CMD to redder colors compared to what is observed in 
the optical would be expected for galaxies with somewhat older average stellar populations 
or with larger reddening due to dust. Similarly, whereas the red sequence is relatively narrow 
in the optical, a color distribution skewed to wards the blue would be expected for early-type 
galaxies with some residual star formation (jYi et al.l 120051 ). The departure of the blue and 
red sequences from a Gaussian would imply that we are beginning to resolve some of these 
effects due to the greater sensitivity of the UV light to both variations in the star formation 
rate and the amount of dust. 

Even though a double Gaussian provides a poor fit to the (NUV — r) .i color distribu- 
tion, the peak of each Gaussian provides a robust estimate of the peak of each sequence in 
each absolute magnitude bin. On the other hand, the width a of each Gaussian should be 
interpreted with caution as it only gives some information about the blue edge of the blue 
sequence and the red edge of the red sequence and does not provide a good representation 
of the entire distribution. The resulting parameters of the Gaussian fits are plotted in Fig- 
ure [151 In the figure the circles and squares give the peaks of the red and blue sequences, 
respectively, whereas the error bars denote the a for each Gaussian. The values of a for the 
red sequence lie in the range 0.3 — 0.5 mag while those for the blue sequence lie in the range 
0.5 — 0.6 mag. Also plotted in Figure [15] are the median photometric errors as a function 
of color for comparison. While the values for a for the blue sequence are significantly larger 
than the photometric errors, the values for the red sequence are comparable to the photo- 
metric errors in the color, indicating that the fall-off of the color distribution on the red edge 
of the red sequence is largely consistent with that expected from the errors. 

We have fit a line to the peak color of the red sequence as a function of absolute 
magnitude with the result (NUV — r) 01 = 1.897 — 0.175M r01 . This fit is plotted as the 
dashed r ed line in Figure M and the solid red line in Figure O T his fit has the sa me slope as 
found by lYi et al.l (120051 ) except with a slightly redder intercept. lYi et al.l (120051 ) analyzed a 
sample of morphologically selected early type galaxies, some of which appear to harbor some 
residual star formation, a fact which would tend to pull their fit somewhat towards the blue 
compared to our color-selected sample. 



Following 



Baldry et al 



(120041 ). we fit the peak color of the blue sequence with the sum 
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of a line plus a tanh function with the result 

(M + 20 32 \ 
r,0 - 1 i8i J . (7) 

This curve is overplotted as the solid blue line in Figure [15j Overall the peak of the blue 
sequence increases in color from (NUV — r) .i = 1.8 at the faint end to ~ 3 for the most 
luminous blue galaxies. The color of the blue sequence changes most rapidly around M ri o.i = 
-20.3. 



3.3. Luminosity Functions as a Function of Color 

3.3.1. M r) o.i luminosity functions 

The galaxy luminosity function varies strongly with color even within each sequence. 
This is illustrated in Figure [161 where the M r) o.i luminosity functions are plotted separately 
for one magnitude wide bins in (NUV — r) .i color ranging from zero to seven. These are 
essentially horizontal cu ts through t he CM P in Figure O We have fit the luminosity function 



in each color bin with a ISchechterl (119761 ) function given by 



$(M) =O.41nlO0*lO- a4 ( M - M *)( a+1 )exp{-lO- a4 ( M - M *)}, (8) 

where <fr*, M*, and a are fit to each sequence by minimizing \ 2 - We determined the la errors 
in each of the parameters using the range of solutions within one of the minimum reduced 
X 2 . The best-fit parameters are listed in Table [6] while M* Q l and a are plotted as a function 
of (NUV - r) .i in Figure EH 

The luminosity function for the very bluest bin with < (NUV — r) .i < 1 is very steep 
with a best-fitting faint end slope of nearly —2, although the errors on M* and a for this bin 
are large and highly correlated. As the color increases through the blue sequence, the faint 
end slope a gradually increases, reaching values of ~ —0.5 in between the two sequences at 
(NUV — r)o.i ~ 4. Although the uncertainties become larger, the slope reaches a slighty 
larger value of ~ for the reddest bin. The value of M* similarly varies systematically with 
color, going from —20.4 in the bluest bins to —20.8 at intermediate colors and finally reaching 



21.1 for the reddest galaxies. Qualitatively similar results were found by iBlanton et al. 



( 120011 ) when determining the r-band luminosity function separated by (g — r) color. 



We have computed the total luminosity density within each color bin from the best-fit 
Schechter parameters as p r ,o.i = Jq°° LQ(L)dL = (p*L*T(a+2). The statistical errors in p were 
determined using the range of values corresponding to those fits within one of the minimum 
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reduced x 2 . For values of the faint end slope a < —2, the integral of the Schechter function 
diverges when integrating all the way down to zero luminosity. With the exception of the 
bluest bin, the luminosity function slopes are significantly larger than —2. However, the la 
confidence interval for a for the bluest bin includes values in the region where a < —2. When 
computing the luminosity density in this bin, we have integrated the luminosity function only 
down to M r) o.i = —12. Since this choice is somewhat arbitrary, the luminosity density of 
the bluest galaxies is more uncertain than for the other color bins. However, we note that 
this luminosity reaches abo ut the limits for whi ch the galaxy luminosity function has been 
determined from the SDSS ( Blanton et al.ll^Ood ). 

The luminosity densities are listed in Table [6] and are plotted as a function of color in 
the top panel of Figure [TS] while the fraction of the total luminosity density within each color 
bin is plotted in the bottom panel of the figure. The largest contribution to the luminosity 
density comes from galaxies with 2 < (NUV — r)o.i < 3 and accounts for pa 30% of the total. 
Galaxies bluer than (NUV — r) .i = 4 together contribute 64% of the luminosity density 
while redder galaxies account for 36%. Adding up the the contribution to the luminosity 
density from each color bin, we obtain a total of logp rj o.i = 26.903 ± 0.030 ergs s _1 Hz -1 
Mpc -3 . This is only slightly larger than the luminosity densit y of log/3 r< n.i = 26.845 ± 0.012 
calculated from a much larger sample of SDSS galaxies by iBlanton et all (l2003bl ). after 
converting to our value for the Hubble constant. 



3.3.2. M NUV fi.i luminosity functions 

Similar to the analysis of the M r .i luminosity functions described in the previous 
section, we have computed M^uvfl.i luminosity functions in one magnitude wide bins of 
(NUV — r)o.i color. As for the r-band, we fit Schechter functions to the distribution within 
each color bin. The luminosity functions are plotted in Figure [19] while the best-fitting 
Schechter function parameters are listed in Tabled The best-fit values for M^ UV01 and a 
are plotted as a function of color in Figure [2171 Qualitatively, the results are similar to that in 
the r-band. The faint end slope a is very steep for the bluest galaxies with a value of —1.8. 
The faint end slope gradually increases with color up to a value of —0.6 at (NUV— r)o.i = 3.5. 
For redder galaxies, the faint end slope remains at about this value until increasing slightly 
again in the reddest bin. The value of M^ uvo l increases dramatically from ~ —20 for the 
bluest galaxies to ~ — 15 for the reddest. 

Similar to the r-band, we have computed the NUV luminosity density Pnuv,o.i within 
each color bin. As before, the luminosity function for the bluest galaxies is near a slope —2 
where its integral diverges. Therefore, for the luminosity density in this bin, we have only 
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integrated the luminosity density down to Mnuv,o.i = — 12. The value of the luminosity 
density as a function of color is plotted in the top panel of Figure [21] while the bottom 
panel shows the fraction of the total luminosity density contributed by the galaxies in each 
color bin. As would be expected, the NUV luminosity density is dominated by the blue 
sequence galaxies. Specifically, « 80% of Pnuv,q.i is coming from galaxies with colors in the 
range 1 < (NUV — r) .i < 3. The bluest galaxies, those with < (NUV — r) .i < 1, only 
contribute ~ 6% to the NUV luminosity density. Galaxies as blue as (NUV — r) 01 ~ are 
difficult to produce using models with smoothly decling star formation histories. Such very 
blue galaxies can be reproduced by models with a star formation burst lasting 1 — 100 Myr, 



with little dust and involving a significant fraction of the mass of the galaxy (jTreyer et al. 



19981 ). Clearly, such dust-free starburst galaxies are relatively rare in the local universe and 
do not contribute much to the total UV luminosity density. The blue sequence could in 
principle include galaxies undergoing large starbursts that have relatively large extinctions. 
However, we argue in §3.5 below that the bulk of the galaxies in the blue sequence do not 
have such extreme dust attenuation. 

Adding up the total luminosity density for galaxies of all colors, we obtain a value 
of log Pnuv,o.i — 25.791 ± 0.029 ergs s" 1 Hz" 1 Mpc~ 3 , a va lue consistent to within the 



uncertainties with the value determined by lWyder et al.1 (120051 ) from an UV-selected sample. 



3.4. Comparison of (NUV — r) 01 with (u — r) .i 
Although there are a number of qualitativ e similarities b etwee n the (NUV — r) .i CMD 



presented here and the (u — r) diagram from iBaldry et al.1 (120041 ). there are a few notable 



differences. As already shown in §3.2, the (NUV — r) .i color distributions are not well fit 
by a double Gaussian function, in contrast to the (u — r) distributions. There is an excess 
of galaxies in between the two sequences above that predicted by the double Gaussian. 
Moreover, the blue sequence appears to merge with the red sequence for the most luminous 
galaxies in the (u — r) diagram while there is still a distinct blue peak visible in Figure [TJ] 
up to M r fi.i = —22.75. In addition, the separation between the blue and red sequences, 
compared to their widths is larger at all luminosities than in the (u — r) distributions. 

The reason for these differences is relatively easy to understand and related to the greater 
sensitivity of the NUV band to changes in the recent star formation rate. To illustrate 
this point, we compare directly the (NUV — r)o.i and (u — r)o.i colors for our sample in 
Figure [22j When generating this figure, we excluded galaxies with w-band photometric 
errors larger than 0.3 mag. For blue galaxies, the two colors are very well-correlated with 
a slope of A(u — r) i/A(NUV — r) .i ~ 0.5. However for galaxies with colors redder than 
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(NUV — r)o.i ~ 3.5, there is a change in slope and the (u — r)o.i color begins to increase 
less quickly with (NUV — r) .i than for bluer (u — r) 01 colors. As a result, galaxies that 
are on the red sequence in the (u — r) .i CMD tend to be more spread out in color in the 
(NUV — r)o.i diagram. 

This behavior is basically that predicted based upon simple galaxy mode ls. In Figure 
[22l we overplot as red circles the locations of a few iBruzual fc Charlotl (120031 ) models with 
an age of 13 Gyr, no dust and solar metallicity. The models are plotted for exponentially 
declining star formation histories with five values of the time constant 7 in the range 0.01 — 7.5 
Gyr -1 . These models are capable of reproducing the locus of data points in Figure 1221 and 
in particular the change in slope for (NUV — r)o.i > 3.5. The main resulting difference in 
the observed CMDs is a greater sensitivity of the UV to small changes in the recent star 
formation rate, especially relevant to galaxies on the red sequence or in between the two 
sequences. For reference, the solid black line in the figure indicates the reddening vector in 
this diagram corresponding to a re ddening in the ionize d gas of E(B — V) gas = 0.5 mag, 
assuming the attenuation law from ICalzetti et all (120001 ) . The reddening vector lies nearly 
parallel to the blue sequence, meaning that dust would tend to simply move galaxies along 
the blue sequence in this diagram. 



3.5. Correcting for Dust 



One of the most important obstacles in interpreting the galaxy CMD is the fact that the 
UV minus optical color of a galaxy is not only affected by the galaxy's star formation history 
but by other physical parameters, primarily the amount of dust and the metallicity. We would 
like to understand how much of the variation in color with luminosity that we observe is due 
to which of these physical parameters. The most reliable method for determining the UV 
attenuation is the far-infrared (FIR) to UV flux ratio because it is almost inde pendent of the 



age o f the stellar population, the dust geometry, or intrinsic dust properties (IGordon et al. 



20001 ). However, the vast majority of our sample have no FIR data available. Thus, we are 
forced to use more indirect methods. We have estimated the effects of dust on the colors of 
galaxies in our sample using two methods, the first usi ng the Balmer lines, and the second 
using the empirical dust-SFH-color relation derived by I Johnson et all (120061 ). 
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3.5.1. Balmer decrement 



We have u s ed th e Hydrogen Balmer line fluxes measured in the SDSS fiber spectra by 
Tremonti et al.l (12004 ) to estimate the attenuation for our sample. The intrinsic ratio of Ha 
to H/3 flux is relativel y independent of t he physical conditions within H II regions and has a 
value of R a /3fi = 2.87 (jOster brockl Il989l ) . For those galaxies with Ha and B./3 emission lines 
detected, we have computed a reddening in the ionized gas E gas (B — V) using 

2.5 log (R a /3/R a p,o) 



E gas {B 



V) 



k{KP) - k{Ra) 



(9) 



where k{H.(3) — A; (Ha) = 1.163 for the Galactic extinction law of ICardelli et al.l (119891 ) 
with Ry = 3.1. We fur ther computed the attenuation in the NUV and r bands using 



the ICalzetti et al.l (120001 ) attentuation law, where A 



NUV 



3.63E gas (B - V) and A r 



1.57E gas (B — V). We have used the line flux uncertainties to compute a statistical un- 
certainty SAnjjv in the resulting attenuation. For those galaxies with no Ha or H/3 lines, 
Anuv < 0, or 5A NUV < 0.5 mag, we have assumed E gas (B — V) — 0. For the remaining 
galaxies, the median NUV attenuation is Anuv — 1-2 mag with 99% of the galaxies lying 
in the range A NUV = — 3 mag. 

There are a few caveats to bear in mind when estimating the attenuation using the 
Balmer lines. The relationship between the R a p measured in the 3" diameter fiber spectra 
and the value for the galaxy as a whole is likely complicated. For many star-forming galaxies, 
the metallicity becomes lower with radius, which we would expect to lead to a decrease in 
the dust attenuation with radius as well. In these cases, the UV attenuation within the 
fiber would be higher than for the galaxy as a whole. On the other hand, some galaxies, for 
example with a strong bulge, have most of their star formation occurring in the outer disk 
of the galaxy. For these cases where there would be weak or no emission lines in the fiber 
spectr a, the integrated UV attenuation would be underestimated. Finally, the attenuation 
law of ICalzetti et al.l (120001 ) was calibrated using observations of starburst galaxies. It has 
been shown that the relation ship between the UV attenuation and UV spectral slope for less 
actiye, mo r e "normal" star-forming galaxies is smaller than predic ted by the starburst results 



(IBellll2002l : Kong et al.ll2004l ; ISeibert et al.ll2005l : Buat et al. 



20051) ■ Whe ther the attenuation 



determined from the Balmer decrement and the ICalzetti et al.l (120001 ) law lead to similar 
overestimates of the attenuation remains unclear. 



The dust-corrected CMD is plotted in Figure [23J The color of the faint end of the 
blue sequence changes little since the reddening for these galaxies is small. The reddening 
increases with luminosity such that the overall trend of color with luminosity disappears 
while leaving the dispersion of the blue sequence relatively unchanged. The dust-corrected 
color of the peak of the blue sequence has a value of (NUV — r) ai « 1.7. 
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3.5.2. The dust- SFH- color relation 



By combining a spectroscopic measure of the star formation hist ory relatively insensitive 
to dust, such as the _D n (4000) index (e.g.. iKauffmann et al.ll2003al ). with UV, optical, and 



FIR fluxes, it is po ssible to separa t e the effects of dust and star formation history on the 



colors of galaxies. IJohnson et al.l (120061 ) have used a sub-sample of SDSS galaxies with 



both UV measurements from GALEX and FIR fluxes from Spitzer to determine the UV 
attenuation, as measured by the FIR to UV ratio, as a function of -D n (4000) and color. 
We have utilized the fits given by IJohnson et al.l (120061 ) to estimate the UV attenuation for 
our entire sample. Specifically, we calculated the FUV atte nuation A FU y as a fu nction of 
-D n (4000) and (NUV — r) .i color from the fits in Table 1 of IJohnson et al.l (120061 ) as 



A FUV = 1.27-1.56{D n (4000)-1.25}+1.35{(JV?7F-r)o.i-2}-1.24{D n (4000)-1.25}{(iV?7F-r)o.i-2)}. 

(10) 

We calculated the NUV and r attenuation us ing equation (ITU|) an d assuming Anuv = 



0.81 Apuv and A r = 0.35A FUV , derived from the lCalzetti et al. 



(120001 ) attenuation curve. 



In order to illustrate the range of NUV attenuations in our sample, in Figure EMI we plot 
the (NUV— r)o.i color as a function of _D n (4000) from the SDSS spectra. The contours in the 
figure reflect the number of galaxies in the sample and have not been weighted by 1/V max . 
Overall, the (NUV — r) .i color and the _D n (4000) index are correlated with redder galaxies 
having on average larger D n (A000) values, indicating older stellar populations, although 
with a relatively large spread. The dashed red lines plotted over the data in Figure EM] are 
lines of constant A^uv using equation (ITU1) . The bulk of the blue sequence galaxies in our 
sample, with absolute magnitudes near M*, have Anuv = 1 — 2 mag, with the mode of 
the distribution at A NUV « 1.3 mag. The attenuation decreases on average with decreasing 
color and reaches values near zero for the bluest galaxies. Since the distribution of galaxy 
colors along the blue sequence do not lie perpendicular to lines of constant attenuation in 
Figure EM] the variation in color along the blue sequence is not entirely due to dust alone — 
the star formation history and metallicity also play a role. 

While the (NUV 

— ?")o.i color is relatively tightly correlated for blue and red sequence 
galaxies, sources with intermediate colors tend to exhbibit a larger spread in D n (4000). 
Indeed, galaxies with these intermediate colors tend to have a large range of attenuation 
ranging from one to three magnitudes. This would indicate that some galaxies are located 
in between the two sequence simply because they are star-forming galaxies with a lot of dust 
reddening while other galaxies are located there due to their older average stellar populations. 



In addition to dust and star formation history, it is important to keep in mind that 
aperture effects can play an important role, particularly for galaxies in between the two 
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sequences because the _D n (4000) measurements only sample the inner 3" of each galaxy. For 
galaxies with strong bulges, for example, the SDSS spectra sample mostly the bulge resulting 
in an old value for D n (4000) whereas the color would tend to indicate a somewhat younger 
age since it measures the flux from both the bulge and any recent star formation in the disk. 
This may explain the population of galaxies in Figure [24] that have values of _D n (4000) that 
would place them on the red sequence but with colors that place them in b etween the red 
and blue sequences. However, as the sample used by iJohnson et al.l (120061 ) is drawn from 
the SDSS, their fits should inherently account for this aperture correction on average. 

The relation in equation fflOl) is an empirical des cription of the correl ations observed 
in the joint SDSS- GALEX-Spitzer sample analyzed by IJohnson et al.l (120061 ). As such, it is 
difficult to interpret directly in terms of physical variables. Presumably, it is the galaxies' 
star formation history which varies along each line of constant UV attenuation shown in 
Figure [24] although the lines also implicitly account for the fact that the Z) n (4000) index 
is measured only w ithin the central regions of each galaxy. In addition, the analysis of 
Johnson et al.l (120061 ) assumes a very simple prescription for converting the FIR/UV ratio to 
an UV attenuation and does not take into account heating due to an older stellar population. 
In such cases, particularly relevant for red sequence galaxies with little or no recent star 
formation, the UV and FIR emission may become de coupled. Indeed, the residuals from the 
fit in equation (fTOl) are larger for red galaxies in the IJohnson et al.l (120061 ) sample. 



Another important caveat to bear in mind when interpreting the lines of constant A^uy 
in Figure [24] is the cuts imp osed to define the sam ple used in determining the fit in equation 
(fit)]) . The sample used by IJohnson et al.l (120061 ) was selected by requiring galaxies to be 
detected both in the UV by GALEX and at 24/im by Spitzer. Thus, their sample will 
not include quiescent early type galaxies with no recent star formation. As can be seen 
in Figure [24] the typical attenuation predicted for red sequence galaxies is A NUV ~ 1.5 
mag. This is most likely an overestimate in many cases since many of these red sequence 
galaxies likely do not contain any residual star formation and are simply red due to their 
older stellar populations. Therefore, care must be taken in interpreting the UV attenuation 
for red sequence galaxies in our sample. 

The (NUV— r)o.i vs. M r;01 diagram corrected for dust using the dust-SFH-color relation 
is shown in Figure [251 In this corrected diagram the ridge line of the blue sequence is 
shifted towards the blue compared to the uncorrected diagram while the width has become 
significantly narrower, particularly at the faint end. In contrast to the Balmer dust corrected 
CMD, a trend of color with absolute magnitude remains. Whereas the ridge line of the 
blue sequence in the uncorrected CMD increases from about (NUV — r)o.i ~ 1.8 to 3 
with increasing luminosity, the blue sequence in the dust corrected diagram increases from 
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(NUV — r)o.i ~ 1.3 to 2.2. Thus, roughly a quarter of the change in color of the blue 
sequence with luminosity is due to dust. The red sequence in the corrected diagram also 
shifts to the blue but has a wider spread in color compared to the uncorrected CMD. As we 
have already argued, it is likely that the attenuation estimated for part of the red sequence 
using equation (TlO|) has been overestimated. In addition, a significant number of galaxies 
with intermediate colors remains. 

Despite th i s cav eat for non-star-forming galaxies, we prefer the relations computed by 
Johnson et al.l (120061 ) as the best measure of the UV attenuation for our sample because 
this method is based upon t he total FIR/UV ra tio, which is the most direct measure of 



the UV light absorbed (e.g. iGordon et al.l |2000| ). Deriving an UV attenuation from the 



Balmer lines requires assuming an attenuation law to convert from the extinction in the 
Balmer lines to that in the UV in addition to the uncert ainty associated with only having 
measurements in the central 3" of each galaxy. Although ICalzetti et al.l (119941 ) showed that 



E s tars(B — V) — QAAE gas (B — V) for starburst galaxies, it remains unclear whether this 
relation pertains to more "normal" galaxies. 

Besides attenuation due to dust and the star formation history, the colors of galaxies can 
also be affected by the metallicity. We have estimated the effect of met allicity on the colors 



of star- forming galaxies using the models of iBruzual Charlotl (120031 ) . For galaxies at an 



age of 12.6 Gyr, no extinction, and exponentially declining star formation histories with time 
constant 7 < 4 Gyr" 1 , the (NUV — r) color varies by 0.2 mag for metallicities between 0.4 
and 2.5 times the solar abundance. For such currently star-forming galaxies, the (NUV — r) 
color becomes bluer with increasing metallicity. This some what counterintuitive result is 



(IBruzual & Chariot 


2003) 


Tremonti et al. 


2004 


). the 



20031 ). Since metallicity tends to increase with galaxy luminosity (e.g., 



decreased somewhat due to this effect. Thus, metallicity variation does not account for the 
remaining change in color along the blue sequence after correcting for dust attenuation. 



3.6. Specific Star Formation Rates as a Function of Stellar Mass 

Although we do not have access to the detailed star formation history of each of the 
galaxies in our sample, we can try to place some basic constraints on their evolution. One of 
the most basic parameters that we may hope to use to constrain the star formation history 
is the specific star formation rate, or the star formation rate per stellar mass SFR/M*. 
We have made estimates of specific star formation rates for our sample as follows. For the 
stellar masses, we have made use of the estimates already included as a part of the MPA/JHU 
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values-added catalogs. As described in detail by lKauffmann et al.l (j2003al ). the stellar masses 
were determined from the z-band luminosities with the mass to light ratio constrained by 
their m o del fi ts to the _D n (4000) and H5a spectroscopic indices. The masses assume a 
Kroupal (120011 ) stellar initial mass function (IMF). We have estimated a star formation rate 
SFR in M yr _1 for eac h galaxy u si ng its dust corrected N UV lumi n osity. We have relied 
upon the relation from iKennicuttl (119981 ) converted to the iKroupal (120011 ) IMF: SFR = 
1.0 x 10~ 28 L v> nuv j where L V>NUV is the NUV luminosity in units of ergs s -1 Hz -1 and SFR 
is the star formation rate in M & yr _1 . 

An alternative method for calculating stellar masses was presented bv lBell et al. d2003h 
Based upon fits of models to galaxies with optical and near-IR photometry, iBell et al.l (120031 ) 



determined relations between the stellar mass-to-light ratio (M*/L) and color for various 



bands. Using their relation betwe en {M*/L) r and 



ing stellar masses with those from iKauffmann et al. 



q — r), we have compared the result- 
(j2003al ). The stellar masses agree on 



ave rage for th e most massive galaxies while those for lower mass galaxies are larger using 
the IBell et al.l (120031 ) relations. The difference in mass increases steadily below 10 10 M Q , 
reaching a difference of ~ 0.4 dex on average a t 10 8 Mm. Throughout th e remainder of this 



discussion, we assume the stellar masses from IKauffmann et al.l ( ]2003al ) but note that the 
specific SFRs for the lowest mass galaxies are likely more uncertain since the origin of this 
difference between the two mass calibrations is not clear. 

The number density of galaxies determined using the V max weighting is shown in Figures 
|2H] and [23 for the two different dust corrections discussed in §3.5. In Figure l2T)j the SFR was 
calculated from the NUV luminosity corrected for dust using the Balmer decr ement, whereas 
the SF R used in Figurel26lwas determined using the dust-SFR-color relation of lJohnson et al. 



(120061 ). The overall trends are similar in both diagrams. The s pecific SFR is correlate d 
strongly with the stellar mass, as previous studies have shown (e.g. iBrinchmann et al.ll2004l ). 
The specific SFR decreases gradually with increasing stellar mass. Above a stellar mass of 
~ 10 10 ' 5 M Q , galaxies with very low specific SFRs begin to dominate, although they are 
present in smaller numbers well below this mass. This is partially a selection effect since our 
sample would be biased against low mass galaxies with low specific SFRs. Indeed, the low 
mass limit in Figures |2T)1 and |2"71 increases as a function of mass due to the UV selection limit 
of our sample. 



There are some notable differences between the specific SFRs in Figures [261 and 1271 The 
width of the blue sequence is substantially broader when using the Balmer decrement dust 
correction. In both diagrams the peak specific SFR in both diagrams is ~ 10~ 10 " 3 yr _1 at the 
high mass end of the blue sequence at 10 11 M . The specific SFR increases somewhat less 
with decreasing mass for the Balmer decrement dust correction than for the iJohnson et al. 
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-9.6 



yr" 



( 120061 ) dust correction. At a stellar mass of 10 8,5 M Q , the peak s pecific SFR is 
for the Balmer decrement dust correction and ~ 10 -9 ' 3 with the iJohnson et al.l (120061 ) dust 
correction. 



There are larger differences in the specific SFRs derived from the two dust corrections for 
red sequence galaxies. The specific SFR is a factor of ~ 10 lower for the Balmer decrement 
dust correction. Determining specific SFRs for red sequence galaxies is problematic in both 
methods. The Balmer line dust correction underestimates the amount of dust absorption in 
cases where most of the star formation is taking place in the outer parts of a galaxy that 
are not sampled by the SDSS fiber spectra. On the other hand, the IJohnson et al.l (120061 ) 
method overestimates the dust absorption in some red sequence galaxies. In particular, dust 
heating by older stellar populations or AGN can contribute to the FIR luminosity to a higher 
degree in red sequen ce galaxies. In addition, the galaxies used to derive the correlations in 
Johnson et al.l (120061 ) were all detected at 24/im, which would bias the sample against truly 
quiescent galaxies with very little dust. Regardless of the dust correction applied, SFRs 
derived from the NUV luminosity for red sequence galax ies are problematic because th e 
NUV band can include light from older stellar populations (lYi et al.ll2005l ; iRich et al.ll2005l ). 
In those cases, the NUV luminosity will overestimate the recent SFR. For these reasons, 
caution must be used when interpreting the specific SFRs for red sequence galaxies. Figures 
1261 and [27] are most useful for investigating the star formation histories of blue seq uence 
galaxies. As we have argued in the previous section, we prefer the IJohnson et al.l ( 20061 ) dust 
correction and thus restrict our discussion to Figure [27J through out the remainder of the 
paper. 

With a few additional assumptions, the specific SFR can be related to the ratio of the 
present to past average SFR as b = SFR/ < SFR >= (SFR/M*)TR, where T is the age 
of the galaxy and R is the fraction of the mass f ormed over the galaxy's lifetime that does 
not eventually get returned to the ISM or IGM (IBrinchmann et al.ll2004l ). A typical value 
of R is ~ 0.5 (IBrinchmann et al.l 12004 ) . Assuming that all galaxies started forming stars 
shortly after the Big Bang, we set T = 13 Gyr. For these choices of T and R, a galaxy with 
6 = 1, i.e. a constant SFR, would have a specific SFR of 10~ 9 ' 8 yr _1 . Thus galaxies with 
specific SFRs larger than this value would have current SFRs above their average while those 
with specific SFRs < 10~ 9 ' 8 yr _1 would have current SFRs below their lifetime average. The 
right hand y-axes of Figures 1261 and 1271 gives the corresponding values of log (6). Under these 
assumptions, blue sequence galaxies with M* < 10 10 M would have star formation histories 
increasing moderately with time, with more massive blue sequence galaxies having formed 
their stars at somewhat declining rate with time. As noted earlier when discussing the CMD 
itself, galaxies undergoing extreme starbursts, i.e. objects with log (6) > 1, are rare in the 
local universe and do not contribute significantly to the SFR density. 
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Quantitatively, the specific star formation ra tes estimated here are somewhat larger 
than those determined by iBrinchmann et al.l (120041 ) from their analysis of the SDSS emission 
lines, especia lly for lower mass galaxies . The origin of this difference is not clear. The SFRs 
estimated by IBrinchmann et al.l (120041 ) were derived largely from the emission lines coming 
from H II regions within the SDSS fiber. Since stars capable of producing an H II region 
are very massive stars with lifetimes up to ~ 10 7 yr, compared to stars with lifetimes up to 
~ 10 8 yr that dominate the UV emission, differences between UV and optical emission line 
based SFRs could be due to short time scale variations in the SFRs of galaxies. This could 
in principle be more important for lower mass galaxies where stochastic fluctuations in the 
SFR would be expected as a fraction of the total mass. In addition to this fundamental 
difference, part of the difference could be due to uncertainties in our method for correcting 
for dust attenuation or to uncertainties in the aperture cor rections applied to the S PSS 
measurements. A more detailed comparison of the S FRs from IBrinchmann et al.l (120041 ) and 
GALEX UV determinations is presented elsewhere (ITreyer et al.l 120071 ). 



Discussion 



Due to the greater sensitivity of the UV minus optical colors to very low specific star 
formation rates as compared to optical colors, our galaxy color-magnitude diagrams have 
revealed a few new features. Similar to optical results, our C MP exhibits prom inent blue 
and red sequences. However, whereas in the 



u ■ 



CMD 



e.g. 



Baldry et al. 



|2004j), the blue 

sequence appears to merge with the red sequence at the luminous end, in the (NUV — r) .i 
CMD, the blue sequence remai ns as a separate peak in the color distr ibution up to iW r o.i — 
-23. Also in the optical CMDs faaldrv et al.lbooi iBalogh et al.ll2004h . the color distribution 
in each absolute magnitude bin is well fit by the sum of two Gaussians. This color bimodality 
naturally leads to an interpretation in which there are simply two populations of galaxies 
with nothing in between. As we have shown in Figure [T2l the color distributions as a function 
Of M r . o.i are not well fit by a double Gaussian function due to an excess of galaxies in between 
the red and blue peaks. These galaxies in the "green valley" between the red and blue peaks 
would suggest that the properties of galaxies are not strictly bimodal and that there really 
exists more of a continuum of properties between star-forming and quiescent galaxies. 

Besides the colors, we have also shown that the r .i-band luminosity function shape 
varies systematically with color. While there is a steady increase in the faint end power law 
exponent a with color across the blue peak, the slope levels off in the "green valley." In 
addition, the decrease in the value of the characteristic magnitude M* 01 appears to level off 
at intermediate colors as well. This change in the luminosity function shape at intermediate 
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colors would also suggest that the " green valley" galaxies coincide with a physical difference 
in the galaxy populations. 

Interpreting the change in the peak color of each sequence with luminosity is complicated 
by the fact that a galaxy's color is affected not only by its star formation history, but also 
its metallicity and dust content. For morphologically selected early type galaxies, there is 
little dust and the change in color is due to some combination of metallicity and age of the 
stellar populations. While the red sequence, as defined by the CMD, contains galaxies with a 
range of morphologies, the dominant effects are still likely the age and metallicity, although 
a significant fr action of early-type galaxies appear to have some residual star formation 



flYi et al J 12005( 1 . 



On the other hand, for blue sequence galaxies, we know that their colors are affected both 
by dust attenuation and their SFH, with metallicity playing a minor role. We have attempted 
to determine the effects of dust in two ways. Usin g the Balmer d e creme nt as observed within 



the SDSS fiber spectra in conjunction with the ICalzetti et al.l (120001 ) attenuation law, we 
found that the trend of the peak color with luminosity disappears, which would suggest that 
the average color along the blue sequence is due mostly to dust with the SFH primarily 
responsible for the dispersion in the color at fixed luminosity. 

Our second method to correct for dust redde ning and attenuation relies upon the empir- 



ical correlations derived by iJohnson et al.l (120061 ). Applying these fits to our entire sample 



we found that only about one quarter of the change in color along the blue sequence is due 
to dust. In this case the dispersion of the blue sequence is small and there is a significant 
variation in the average color with luminosity, indicating a decrease in the average galaxy 
age with decreasing luminosity. Due to the requirement that the galaxies used to calibrate 
these relations must be detected in both the UV and at 24//m, there is an inherent bias 
against truly quiescent galaxies with little star formation and dust. Thus, the attenuation 
for the red galaxies is likely overestimated. Nevertheless, we argued in §3.5 why we prefer 
this dust correction. 

While the detailed morphology of the CMD depends upon which dust correction is 
assumed, there remain galaxies in the "green valley" between the two sequences in both 
versions of the dust-corrected CMD. While there are indeed some galaxies with intermediate 
colors that are simply very reddened versions of blue sequence galaxies, many do have SFHs 
weighted to older ages. We would expect that as time progresses, the red sequence is building 
up as galaxies exhaust their gas and cease forming stars. With the data presented here we 
may be beginning to separate out the population undergoing this transition. Of course, it is 
important to note that just based upon the integrated measurements here it is not possible 
to infer the evolution of different structures within each galaxy. In fact, galaxies with bulges 
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would be expected to lie at intermediate colors, where the SFH of the bulge may be quite 
different than that of the disk. In these cases, the integrated color would reflect the relative 
strengths of the bulge and disk. In addition, it is not possible just based upon the simple 
diagnostics employed in this paper to infer whether galaxies in the green valley are star- 
forming galaxies transitioning to the red sequence for the first time or whether they were 
already on the red sequence and underwent a burst of star formation after merging with a 
gas rich galaxy. The rate at which galaxies ar e moving from the blue to the red sequence is 



explored in more detail in iMartin et al.l (120071 ). 



A remarkable feature of the specific SFR as a function of stellar mass in Figure [271 is the 
relatively small spread of a factor of 2 — 3 in SFR/M* at a given stellar mass. It is difficult 
for galaxies to have a large SFR/M* ratio for an extended period of time due to exhausting 
the available gas supply as well as feedback from supernovae which tends to heat the gas 
and prevent it from forming stars. On the other hand, the energy input from supernovae 
can help maintain lower levels of star formation over an extended period of time as long as 
the halo is massive enough to retain the heated gas. 

The behav i or see n in Figure [57J is broadly consistent with the theoretical models of 



Cattaneo et al.l (120061 ) which succeeded in reproducing the distribution of galaxies in the 
optical (u — r) CMD. In their models, galaxies below a critical halo mass of 10 12 M accrete 
gas from cold streams. This gas fuels star formation that is self-regulated due to feedback 
from supernovae and stellar winds. These lower mass galaxies gradually gain in mass, move 
up in luminosity along the blue sequence, and become redder. The gradual decrease in the 
specific star formation rate with stellar mass seen in Figure [27J would be consistent with the 
predictions of this model. The gradual transition to the dominance of red galaxies above a 
stellar mass of ~ 10 10 ' 5 M Q would correspond in this picture to the critical halo mass of 10 12 
M Q above which star formation is quenched due to shock heating of the gas and feedback 
from AGN. In their model lower mass galaxies can also have their star formation quenched 
if they happen to become satellites within a halo above the critical mass. This tends to 
populate galaxies along the fainter end of the red sequence and serves to broaden somewhat 
the mass range over which the transition from star-forming to quiescent galaxies occurs. The 
most luminous blue galaxies help constrain the value of the critical halo ma ss M s hock above 



which gas is not able to cool and form stars. While ICattaneo et al.l (120061 ) have used the 
(u — r) CMD to constrain the value of M shock, our observations of blue sequence galaxies at 
slightly larger luminosities may correspond to somewhat larger values for M shock- 



In their analysis of the evolution of galaxies in the (U—B) CMD out to z = l. lFaber et al. 



( 120061 ) discussed the origin of red sequence galaxies and how blue sequence galaxies transition 



from blue to red. In their favored scenario, galaxies first grow moderately in mass along the 
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blue sequence, then undergo a merger which halts further star formation, and finally undergo 
a series of gas-free mergers that move the galaxies up the red sequence. Such a scenario is 
broadly consistent with the results presented here. The excess of galaxies in between the two 
sequences would in this case be consistent with galaxies in the process of turning off their 
star formation and making the transition to the red sequence. As can be seen in Figures 
[T6l and [TTJ, the M rj0 .i luminosity functions for the bluest galaxies have a steeper slope and 
fainter M* than for galaxies in between the two sequences. If large numbers of lower mass 
galaxies were stopping star formation and moving towards the red sequence, we would expect 
to see a steeper faint end slope to the luminosity function for galaxies with intermediate and 
red colors. Thus, the luminosity functions are more consistent with a scenario in which the 
ancestors of galaxies on the red sequence are weighted towards the luminous end of the blue 
sequence. The bright end cut-off of the r-band luminosity functions is remarkably similar 
for galaxies with 2 < (NUV — r) .i < 5 while there appears to be a relatively sharp jump 
in the bright end at (NUV — r)o.i ~ 5. Thus, the most massive galaxies are unlikely to be 
the descendants of the merging of galaxies at the luminous end of the current blue sequence. 
Either these galaxies are the descendants of much more massive blue sequence galaxies that 
are not present in the nearby universe, or they are the result of "dry" mergers of smaller 
mass red sequence galaxies, as iFaber et al.l (120061 ) argued. 



5. Summary 

We have determined the volume density of galaxies in the local universe as a function 
of absolute magnitude M r01 and (NUV — r) .i and (FUV — r) .i colors based upon a 
sample of galaxies observed in the UV by GALEX and with optical data from the SDSS. 
The galaxies in these CMDs separate into well-defined blue and red sequences that become 
redder with increasing luminosity. While the most luminous galaxies are on the red sequence, 
a separate blue peak is detectable as bri ght as M r< n.i ~ ~23 . In contrast to CMDs relying 



solely on an optical color such as (u — r) (IBaldry et al.ll2004l ). the color distribution at fixed 



absolute magnitude is not well fit by the sum of two Gaussians due to an excess of objects at 
intermediate colors between the blue and red peaks. The greater separation between the blue 
and red sequences is a consequence of the greater sensitivity of the UV bands to very low 
levels of recent star formation. The ro.i-band luminosity function shape varies systematically 
with color with the faint end slope a gradually increasing across the blue sequence, reaching 
a value of a ~ —0.6 at intermediate colors before increasing even more for the reddest 
galaxies. We have used these fits to the luminosity functions to derive the fraction of the 
luminosity density in the local universe as a function of color. Dust-free starburst galaxies 
with colors (NUV — r)o.i < 1 are rare in the local universe and account for only about 5% 
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of the NUVq.1 luminosity density. About 80% of the NUVq.i luminosity denisty is emitted 
by blue sequence galaxies with colors 1 < (NUV — r) 01 < 3. 



We have used both the Balmer decrement and the dust-SFH-color relation of Johnson et al. 



(120061 ) to estimate the effect of dust on the galaxy colors and absolute magnitudes. For 
the Balmer decrement method, the increase in color with luminosity along the blue se- 
quence is due entirely to dust with the dispersion at fixed absolute magnitude relatively 
unc hanged. On th e other hand, the blue sequence color in the CMD corrected for dust using 
the iJohnson et al.l (120061 ) method does still increase with luminosity, indicating that part of 
this change in color is due to t he sta r formation history and not to dust alone. We argue that 
we prefer the IJohnson et al.l (120061 ) method as it is ultimately based upon an attenuation 
derived from the FIR/UV ratio. Regardless of which dust correction we employ, however, a 
significant number of galaxies remain at colors in between the two sequences, indicating that 
not all of the galaxies there are simply dusty versions of blue sequence star-forming galaxies. 



We have used the NUVq.i luminosities corrected for dust u sing thelJohnson et al.l (120061) 
method in conjunction with the stellar masses determined by iKauffmann et al.l (j2003al ) to 
plot the density of galaxies as a function of specific star formation rate SFR/M* and stellar 
mass M*. The dispersion in SFR/M* is only a factor of 2 — 3 at a fixed stellar mass along the 
blue sequence. The value of SFR/M* decreases from values of ~ 10~ 9 ' 3 yr _1 at Af* = 10 8 ' 5 
M to « 10~ 10, 3 near the tip of the blu e sequence at M* = M 11 M . Similar to previous 
optical results ( Kauffmann et al. 2003al lbl). galaxies with low specific star formation rates 
begin to dominate above a stellar mass of about 10 10,5 M Q . 

In addition to the small number of galaxy properties explored here, many galaxies in 
our sample con tain many other me asurements mainly from the SDSS. In a companion paper 
in this volume, iMartin et al.l (120071 ) have estimated the mass flux of galaxies from the blue to 
the red sequence and have discussed some of t he other properties o f the g alaxies in between 
the red and blue sequences. In another paper, ISchiminovich et al.l (120071 ) have investigated 
the correlation of morphology and other characteristics with position in the CMD. While 
detecting red sequence galaxies out to significant distances in the rest-frame UV is very 
difficult, it should be possible to use data from GALEX deep exposures in conjunction with 
ground-based photometry and spectra to investigate the evolution of the blue sequence with 
redshift. In addition, the variation of the CMD with local galaxy density should provide 
interesting contraints on the nature of the galaxies in between the blue and red peaks as well 
as models of the physical processes affecting the evolution of galaxies in the CMD. 



GALEX (Galaxy Evolution Explorer) is a NASA Small Explorer, launched in April 
2003. We gratefully acknowledge NASA's support for construction, operation, and science 
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analysis for the GALEX mission, developed in cooperation with the Centre National d'Etudes 
Spatiales of France and the Korean Ministry of Science and Technology. 

Facilities: GALEX 
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Table 1. Sample selection limits 



Parameter 



Limits 



SDSS primary target flag 

r-band magnitude 3, 
r-band magnitude error 
redshift 
redshift confidence 
GALEX field radius 
GALEX exposure time 
NUV flag 
FUV and NUV apparent magnitude 3 



TARGET_GALAXY 
14.5 < r < 17.6 

O r < 0.2 

0.01 < z < 0.25 

Zconf > 0.67 

fovjradius < 0?55 

* > 1000 s 
nuv-artifact < 1 
16.0 < N(F)UV < 23.0 



a The r-band magnitude limits were applied after correcting 
for foreground Galactic extinction while those for the NUV and 
FUV were applied before the extinction correction. 
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Table 2. Galaxy volume densities vs. (NUV — r)o.i and M r> o.i a 



a The data for this table are available in the electronic edition of the journal only. The 
data are in units of Mpc -3 mag -2 . 



Table 3. Errors in the galaxy volume densities vs. (NUV — r) .i and M r ^ i 



a The data for this table are available in the electronic edition of the journal only. The 
data are in units of Mpc -3 mag -2 and correspond to the statistical errors only calculated 
using equation (Q. 



Table 4. Galaxy volume densities vs. (FUV — r) .i and M r o.i a 



a The data for this table are available in the electronic edition of the journal only. The 
data are in units of Mpc -3 mag -2 . 



Table 5. Errors in the galaxy volume densities vs. (FUV — r)o.i and M ri o.i 



a The data for this table are available in the electronic edition of the journal only. The 
data are in units of Mpc -3 mag -2 and correspond to the statistical errors only calculated 
using equation (jBj). 
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Table 6. r .i Schechter Function Parameters vs. Color 



(NUV — r )o.i 


log0* 


M rfiA 


a 


log Prfi.l 




(Mpc- 3 ) 






(ergs s _1 Hz -1 Mpc -3 ) 


0.50 


-4.603 ± 0.960 


-20.454 ±0.995 


-1.913 ±0.840 


24.931 ±0.838 a 


1.50 


-2.871 ±0.185 


-20.331 ±0.250 


-1.465 ±0.175 


26.121 ±0.084 


2.50 


-2.564 ± 0.060 


-20.786 ±0.105 


-0.894 ±0.140 


26.368 ±0.032 


3.50 


-2.775 ± 0.045 


-20.711 ±0.125 


-0.357 ±0.220 


26.103 ±0.035 


4.50 


-2.962 ±0.065 


-20.874 ±0.155 


-0.579 ± 0.235 


25.975 ± 0.042 


5.50 


-2.864 ± 0.070 


-21.158 ±0.155 


-0.596 ±0.185 


26.187 ±0.039 


6.50 


-3.474 ±0.125 


-21.127 ±0.345 


+0.022 ±0.660 


25.621 ±0.087 


All colors 








26.898 ± 0.022 b 



a When determining p^o.i for the bluest bin, we have only integrated the luminosity function 
down to an absolute magnitude M r)0 .i = —12, rather then to zero luminosity in the other color 
bins. 

b The total luminosity density is derived by summing all of the values for p r ,o.i in each color 
bin. 

Note. - All the errors quoted in this table are the statistical errors only derived by taking 
all combinations of parameters within one of the minimum reduced x 2 ■ All values quoted in this 
table assume H = 70 km s" 1 Mpc~ 3 , fl m = 0.3, and = 0.7. 
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Table 7. NUV .i Schechter Function Parameters vs. Color 



(NUV-r) 0A 


log0* 


1V1 NUV,0.1 


a 


logPAW,0.1 




(Mpc- 3 ) 






(ergs s" 1 Hz -1 Mpc -3 ) 


0.50 


-4.604 ± 0.900 


-19.792 ±0.995 


-1.786 ±0.730 


24.467 ±0.402 a 


1.50 


-2.707 ±0.130 


-18.428 ±0.185 


-1.284 ±0.130 


25.409 ±0.047 


2.50 


-2.566 ±0.050 


-18.341 ±0.090 


-0.943 ±0.105 


25.397 ±0.024 


3.50 


-2.862 ± 0.045 


-17.545 ±0.105 


-0.599 ±0.145 


24.744 ±0.026 


4.50 


-3.047 ±0.060 


-16.605 ±0.125 


-0.746 ±0.150 


24.192 ±0.030 


5.50 


-2.880 ± 0.045 


-15.689 ±0.095 


-0.637 ±0.120 


23.985 ±0.026 


6.50 


-3.469 ± 0.080 


-15.039 ±0.255 


-0.171 ±0.405 


23.159 ±0.055 


All colors 








25.791 ±0.029 b 



a When determining Pnuv,o.i for the bluest bin, we have only integrated the luminosity function 
down to an absolute magnitude Mnuv&a = — 12, rather then to zero luminosity in the other color 
bins. 

b The total luminosity density is derived by summing all of the values for Pnuv,o.i in each color 
bin. 

Note. - All the errors quoted in this table are the statistical errors only derived by taking 
all combinations of parameters within one of the minimum reduced x 2 - All values quoted in this 
table assume H = 70 km s" 1 Mpc" 3 , VL m = 0.3, and = 0.7. 
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0.00 0.05 0.10 0.15 0.20 0.25 0.30 

Redshift 

Fig. 1. — The redshift distributions for the NUV sample. The solid black line indicates the 
redshift distribution for the SDSS main galaxy sample covered by the GALEX fields with 
NUV exposure times greater than 1000 sec and satisfying the SDSS selection criteria in 
Table HJ The dashed red line is the redshift distribution for galaxies with a GALEX match 
within 4" and satisfying the GALEX plus SDSS selection criteria in Table [U 
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Fig. 2. — The same as in Figure [T] except for the FUV sample. 
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Fig. 3. — The apparent r-band magnitude as a function of a) (NUV — r) color for the 
7V[/V-selected sample and b) (FUV — r) color for the Ff/V-selected sample. In both panels 
the horizontal dashed lines indicate the r-band bright and faint magnitude limits while the 
diagonal dashed lines indicate the FUV and NUV magnitude limit of 23, corrected to the 
median Galactic extinction in the samples. In the NUV sample, the lack of galaxies in the 
upper right hand corner would indicate that the fall-off in galaxies at the red end of the color 
distribution is real and not due to the UV magnitude limit. This does not appear to be the 
case for the FUV sample, where the data go all the way up to the magnitude limit. In this 
case, we may not be sampling the full color distribution of galaxies. 
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Fig. 4. — The fraction of galaxies in the SDSS main galaxy sample with GALEX NUV (red 
solid line) and FUV (blue dashed line) matches as a function of r magnitude. The fall in the 
NUV completeness for r < 16.5 is due to the loss of red galaxies below the NUV magnitude 
limit. The completeness at the bright end is about 90%. While most of these galaxies are 
visible in the GALEX images, they are not matched to the SDSS sample due to shredding 
of large galaxies in the GALEX data or due to the UV centroid being offset from the SDSS 
position by more than the 4" search radius used to match the two catalogs. The fall-off in 
the FUV completeness occurs at r > 15.5 due to the larger spread in {FUV — r) colors 
compared to the (NUV — r) distribution. 
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Fig. 5. — The fraction of galaxies in the SDSS main galaxy sample with GALEX NUV 
matches as a function of (g — r) 01 and M r)01 . The completeness for optical blue sequence 
galaxies is greater than 90% while the completeness then drops for galaxies with colors 
id — r )o.i > 0.8. For galaxies along the red sequence, the completeness in the NUV is 
between 30 — 50%. This incompleteness is due to the UV magnitude limit and is accounted 
for when calculating the volume densities shown in the figures below. 
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Fig. 6. — Similar to Figure [5] except for the FUV sample. The fraction of SDSS main galaxy 
sample galaxies with GALEX FUV matches is about 90% for blue galaxies and drops to 
values of 10 — 30% along the red sequence. 
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Fig. 7. — The distribution of galaxies in the (NUV — r)o.i vs. M rj o.i plane. The contours 
are spaced logarithmically from a minimum of five galaxies per bin up to 300 per bin. Below 
the minimum contour level, the positions of individual galaxies are plotted. The median 
photometric errors as a function of color are shown along the left-hand side of the plot. 
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Fig. 8. — The same as Figure [7] except for (FUV — r) .i vs. M r)0 .i. The data are plotted in 
bins that are 0.2 mag wide in color and 0.5 mag wide in absolute magnitude. The contours 
are spaced logarithmically from a minimum of five galaxies per bin up to a maximum of 250 
per bin. The median photometric errors as a function of color are shown along the left hand 
side of the plot. 
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Fig. 9. — The volume density of galaxies as a function of (NUV — r) .i and M r)0 .i using the 
V max method. The density is determined in 0.2 mag wide bins in color and 0.5 mag wide 
bins in absolute magnitude. There are twelve logarithmically spaced contours ranging from 
10~ 5 ' 5 to 10~ 2 ' 3 Mpc -3 mag -2 . The dashed lines indicate the ridge lines of the red and blue 
sequences as determined from the Gaussian fits shown in Figure [15] below. 
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Fig. 10. — The volume density of galaxies as a function of (FUV — r)o.i and M r> o.i using the 
Vmaz method. The density is determined in 0.2 mag wide bins in color and 0.5 mag wide 
bins in absolute magnitude. The contours are logarithmically spaced from 10~ 5 ' 5 to 10~ 2 ' 3 
Mpc -3 mag~ 2 . 
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Fig. 11. — The volume density of galaxies as a function of (NUV — r) .i and Mnuv,o.i using 
the V max method. The density is determined in 0.2 mag wide bins in color and 0.5 mag 
wide bins in absolute magnitude. The contours are logarithmically spaced from 10~ 5 5 to 
10 -2 ' 3 Mpc -3 mag -2 . This is the same sample as in Figure [9] except plotted as a function 
of the NUV absolute magnitude. The sample reaches significantly fainter NUV absolute 
magnitudes for the red galaxies due to the SDSS r-band magnitude limit. 
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Fig. 12. — The volume density of galaxies from the V max method as a function of (NUV — 
r) .i color in 0.5 magnitude wide absolute magnitude bins. The midpoint of each absolute 
magnitude bin is labeled in the upper right of each figure. The blue and red dashed lines 
mark the limits used when fitting Gaussians to the red and blue sequences separately. The 
solid black line is the sum of the two Gaussians. Note that in most absolute magnitude bins, 
the color distribution is not well fit by the sum of two Gaussians as there is a significant 
excess of objects in between the two sequences. 
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Fig. 13. — The continuation of the previous figure. 
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Fig. 14. — The continuation of the previous figure. 
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Fig. 15. — Plot of the Gaussian fits to the red and blue sequences in Figures [T2l— TT4l The red 
circles and blue squares refer to the means of the Gaussians for the red and blue sequences, 
respectively, while the error bars on each point refer to the width a for the fit to each 
sequence. The solid red line is a simple linear fit to the red sequence while the dashed red 
line indicates the fit to the color as a function of absolute magnitude from lYi et al.l (120051 ) 
for morphologically selected elliptical galaxies. The solid blue line is the best-fit T-function 
to the blue sequence. The median photometric errors as a function of color are shown along 
the left hand side of the plot. 
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Fig. 16. — The M r) o.i luminosity functions separated in one magnitude wide bins of (NUV — 
r)o.i color. The errors bars represent the statistical error derived from the V max calculation. 
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Fig. 17. — The best-fit values for M* 01 (top panel) and a (bottom panel) as a function of 
(NUV — r)o.i color from Schechter function fits to the luminosity functions plotted in Figure 
[T6l The error bars are derived from the range of parameters within one of the minimum 
reduced x 2 . 
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Fig. 18. — (a) The r .i band luminosity density p r ,o.i as a function of (NUV — r) .i color 
calculated by integrating the luminosity functions using the best-fit Schechter function pa- 
rameters listed in Table [6j (b) The fraction of the total luminosity density contributed from 
each color bin. 
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Fig. 19. — The M NUVt o.i luminosity functions for galaxies in one magnitude wide bins in 
color derived using the V max method. 
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Fig. 20. — The best-fit values for M^ UV01 (top panel) and a (bottom panel) as a function of 
(NUV — r)o.i color from Schechter function fits to the luminosity functions plotted in Figure 
[T9l The error bars are derived from the range of parameters within one of the best-fit y 2 . 
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Fig. 21. — (a) The NUVq.\ band luminosity density Pnuv,q.i as a function of (NUV — r)o.i 
color calculated by integrating the luminosity functions using the best-fit Schechter function 
parameters listed in Table [7J (b) The fraction of the total luminosity density contributed 
from each color bin. 
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Fig. 22. — The number of galaxies (not weighted by l/V max ) as a function (u — r) 01 and 
(NUV — r)o.i- These are the same galaxies plotted in Figure [7] except that galaxies with 
w-band errors greater than 0.3 mag were excluded. The data are plotted as contours where 
the density of galaxies is high and as individual points where it is low. The solid red circles 



arc 



Bruzual fc Charlotl (120031 ) models at an age of 13 Gyr with no dust, solar metallicity, 



and exponentially declining star formation histories. The time constant in units of Gyr -1 in 
the star formation history for each model is indicated at the bottom of the figure. The black 
line indicates t he red dening vector in this diagram for E(B — V) gas = 0.5 and assuming the 
Calzetti et al.l (120001 ) attenuation law. 
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Fig. 23. — The number density of galaxies as a function of dust-corrected color (NUV — 
r )o.i,corr and absolute magnitude M n o.i. The NUV and r band attenu ation has been calcu- 
lated from the Balmer decrement and assuming the extinction law of ICalzetti et al.l (I1994I ) 
as described in §3.5.1. The contours are spaced logarithmically from 10 -5 ' 5 to 10~ 2 ' 3 Mpc -3 
mag -2 . 
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Fig. 24. — The number of galaxies (not weighted by 1/V max ) as a function of (NUV — r) 01 
and _D n (4000) from the SDSS spectra. The data are plotted as contours where the density 
of galaxies is high and as individual points where it is low. The dashed curves indicate lines 
of constant NUV ex tinction ranging from 0.2 to 3 magnitudes as derived from the fits in 



Johnson et all (12009 ). 
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Fig. 25. — The volume density of galaxies as a function of dust-corrrected (NUV — r) .i,corr 
and M r) o.i jCorr . The FUV attenuati o n was derived using equation QUI) in the text based 



upon the analysis of Ijohnson et al.l (120061 ). The corresponding val ues of Anjjv and A T 



were calculated using Anuv = 0.81Apuv and A r = 0.35Apuv from the lCalzetti et al.l (120001 ) 
attenuation law. The contours are spaced logarithmically from 10~ 5,5 to 10~ 2 ' Mpc" 
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Fig. 26. — The volume density of galaxies as a function of specific SFR and stellar mass 
M*. The SFR has been calculated from the NUV luminosity, corrected for dust using the 
Balmer lines measured from the SDSS fiber spectra as described in the text. The density 
was calculated in bins 0.2 dex wide in mass and 0.1 dex wide in specific star formation rate. 
The contours are spaced logarithmically from 10~ 5 to 10~ L9 Mpc~ 3 dex~ 2 . We con verted 
the NUV luminosities to SFRs using the conve rsion fac t or giv en in iKennicuttl (119981 ) after 
applying a correction factor to convert to the iKroupal (120011 ) IMF used to calulcute the 
stellar masses. The axis on the right hand side gives the value of the logarithm of the ratio 
of current to past averaged SFR, log (6), calculated as described in the text. 



-64- 



-8 




7 8 9 10 11 12 

Log M. (IVy 



Fig. 27. — Similar to Figure [26], except that the SFR has been calcula ted from the NUV 
luminosity using the dust-SFH-color relation from I Johnson et al.l (I2006I ). The contours are 
spaced logarithmically from 10~ 5 to 10 -1 ' 5 Mpc -3 dex -2 . 



